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Elliptic  Length  Scales  in  Laminar,  Two-Dimensional 

Supersonic  Flows 


James  H.  Miller1 

Air  Force  Research  Laboratory,  Wright-Patterson  AFB  45433 

Nomenclature 

Cp  =  pressure  coefficient,  Cp  =  1  p  Pc°  2 

y  / 2  Poo  ^00 

/  =  length  of  elliptic  region,  m 

L  =  reference  length,  m 

j-  =  ratio  of  upstream  influence  length  to  undisturbed  (flat  plate)  boundary  layer  thickness 

j-  =  ratio  of  elliptic  region  length  to  undisturbed  (flat  plate)  boundary  layer  thickness 

/„  =  upstream  influence  length,  m 

ld  =  downstream  influence  length,  m 

V 

M  =  Mach  number,  ratio  of  velocity  to  speed  of  sound,  M  —  , 

V Ip 

p  =  pressure,  N/m2 

Re  =  Reynolds  number,  ratio  of  inertia  forces  to  viscous  forces.  Re  —  Po°v°°  L 

lloo 

t  =  time,  ^ 

u,v,w  =  Cartesian  velocity  components  in  the  x,y,z  directions  respectively,  m/s 

Vx  =  reference  or  freestream  velocity,  m/s 

x,y,z  =  Cartesian  coordinates,  m 

Ap  =  reference  change  in  pressure  (assumed  positive)  relative  to  freestream  (e.g.,  pc  -  px%  where  pc  is  the 
inviscid  pressure  downstream  of  a  shock  interaction  with  a  flat  surface  or  compression  ramp),  N/m2 
SL  =  boundary  layer  thickness  at  location  L,  m 

y  =  ratio  of  specific  heats,  1 .4  for  air. 

9  =  compression  ramp  angle,  degrees 

1  Lead  Senior  Aerospace  Engineer,  Aerospace  Systems  Directorate,  AFRL/RQHV,  B18A,  Rm  A114,  1950  5th 
Street,  Wright-Patterson  AFB  OH,  45433. 
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t  =  shear  stress,  N/m2 

ju  =  viscosity,  kg/(m  s) 

Superscripts: 

~  =  nondimensional  quantity  normalized  by  freestream  conditions  or  reference  change  in  pressure 

Subscripts: 

d  =  downstream  boundary  of  elliptic  region 

L  =  x  location  of  compression  ramp  leading  edge  or  shock  impingement  point 
u  =  upstream  boundary  of  elliptic  region 

xx, xy  =  Cartesian  components  of  shear  stress 
oo  =  freestream  conditions 


I.  Introduction 

Several  theories  [1-3]  and  correlations  [4]  exist  to  describe  the  length  scales  associated  with  the  interaction  of  a 
pressure  disturbance  with  a  boundary  layer  in  a  predominantly  supersonic  flowfield.  The  theoretical  approaches  are 
essential  to  understanding  the  detailed  structure  of  the  viscous  layer  as  it  encounters  pressure  disturbances  that  are 
either  compressive  or  expansive  in  nature.  The  interactions  in  compressive  flows  are  more  straightforward  for 
laminar  flows  but  take  on  different  characteristics  when  applied  to  turbulent  interactions,  requiring  that  theoretical 
descriptions  become  more  sophisticated.  Throughout  the  literature,  the  length  scales  of  these  interactions  are 
typically  described  in  terms  of  the  undisturbed  boundary-layer  thickness;  and,  the  correlations  developed  to  describe 
these  length  scales  take  many  forms  based  on  dimensional  analysis  and  geometric  considerations.  In  the  present 
article,  another  approach  to  these  types  of  theories  and  correlations  based  on  a  similarity  perspective  is  presented. 
New  correlations  for  laminar  compressive  interactions  for  two-dimensional  compression  ramps  and  shock 
impingements  for  adiabatic  walls  have  been  developed.  These  correlations  highlight  a  more  universal  form  which  is 
based  on  a  new  similarity  parameter  of  the  Navier-Stokes  equations.  Results  show  that  compressive  interactions  can 
be  combined  into  a  single  correlation  for  both  upstream  influence  lengths  and  total  elliptic  lengths  of  interactions. 

II.  Analysis 

Pressure  and  viscous  forces  play  a  dominant  role  in  determining  the  length  scale  of  an  elliptic  region  caused  by  a 

pressure  disturbance  in  viscous  supersonic  flow.  Elliptic  regions  are  defined  as  regions  in  which  mathematically 
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elliptic  methods  (time-iterative)  must  be  used  to  accurately  capture  the  details  of  the  flowfield,  while  parabolic 


regions  can  be  solved  using  space-marching  methods  where  no  significant  upstream  influence  affects  the  flow.  The 
ratio  of  inertial  forces  to  viscous  forces  in  the  form  of  Reynolds  number  plays  an  important  role  in  defining  the 
length  scales  of  the  elliptic  regions  because  of  the  importance  of  the  boundary  layer  thickness.  In  this  work  another 
similarity  parameter  is  utilized  to  determine  the  length  of  elliptic  regions,  and  this  similarity  parameter  represents  the 
ratio  of  pressure  forces  to  viscous  forces.  To  the  author’s  knowledge,  there  are  presently  no  similarity  parameters  for 
the  Navier-Stokes  equations  that  have  been  described  in  the  literature  for  shock  boundary  layer  interactions  where 
the  ratio  of  pressure  forces  to  viscous  forces  is  an  important  consideration.  Although  reference  to  a  similar  parameter 
was  discussed  in  Ref.  [5]  to  capture  small  scale  disturbances,  the  present  work  is  focused  on  the  larger  scale  of  the 
interaction. 

This  discussion  is  restricted  to  the  2D  compressible  momentum  equation  in  the  Cartesian  x  direction[6] : 
dpu  d  d 

~dT  +  dx^u‘ +  p  ~  Txx]  +  dp vuv  ~  Xxy'  =  ° 

with 


^ XX 


dv  \  f  du  dv  \ 

)  ’  Txy  =  +  dx) 


In  the  following  nondimensionalization  approach,  nondimensional  quantities  are  annotated  with  a  tilde: 


x,  y 


x  y 

L’L  ’  L 


'/l4 


;  p 


;  u ,  v 


u  v 


;  p 


P_ 
A  p 


There  is  nothing  unusual  about  this  approach  to  nondimensionalization,  except  for  the  pressure.  In  this  approach  the 
pressure  is  nondimensionalized  by  a  change  in  pressure  that  is  representative  of  the  interaction  causing  the  formation 
of  the  elliptic  region.  For  example,  an  appropriate  value  for  A p  would  be  the  overall  pressure  rise  associated  with  a 
2D  compression  ramp  in  an  inviscid  flow.  For  internal  flows,  the  pressure  change  could  be  tied  to  the  difference 
between  a  nozzle  exit  pressure  and  a  back  pressure.  For  rather  simple  2D  external  flows,  the  change  in  pressure  is 
easily  computed  a  priori  using  analytical  shock  relations;  however,  for  3D  interactions,  the  change  in  pressure 
would  likely  need  to  be  computed  with  more  sophisticated  computational  fluid  dynamics  (CFD)  methods. 
Additionally,  for  3D  interactions,  the  length  scales  would  require  determination  in  spanwise  as  well  as  streamwise 
directions. 
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Substituting  for  the  dimensional  quantities  in  the  momentum  equation,  and  then  rearranging  terms  and  simplifying 
leads  to  the  following  equation  in  nondimensional  form: 


dpu  d  ,  1  /  Re  Ap  \ 

-W  +  +  TeivMiz)  p  -  11  + 


—  \puv  -  —fxy\  =  0 


In  preliminary  studies  of  scaling  elliptic  regions,  it  is  more  useful  to  treat  the  parameter  multiplying  the  pressure 
term,  R^Ap  ,  as  independent  of  Reynolds  number  when  there  is  an  elliptic  interaction.  This  approach  is  consistent 

with  the  assumptions  in  boundary-layer  theory  in  which  the  pressure  distribution  is  imposed  on  the  boundary  layer. 
However,  the  present  approach  allows  for  more  insight  into  the  relative  importance  of  the  pressure  term  when 
compared  to  the  viscous  terms  in  the  equation.  From  a  mathematical  perspective,  the  eigenvalue  analysis  of 
Vigneron  et  al.[7]  resulted  in  a  relation  for  the  mathematical  character  of  the  parabolized  Navier-Stokes  equations 
that  utilized  a  weighting  factor  on  the  streamwise  pressure  gradient.  In  that  analysis  of  the  local  mathematical 
character,  the  weighting  factor  was  determined  to  be  a  function  of  the  streamwise  Mach  number  only.  Discussed 
herein  is  how  a  new  similarity  parameter  can  be  used  to  predict  length  scales  of  elliptic  interactions,  but  not  used  to 
explore  details  of  the  local  mathematical  character  of  the  equations  nor  for  an  order-of-magnitude  analysis  of  the 
terms  in  the  equation. 

Before  a  discussion  of  the  correlation  results  can  be  fully  appreciated,  a  description  of  how  the  new  similarity 
parameter  represents  a  ratio  of  pressure  forces  to  viscous  forces  is  warranted: 


Considering  a  fluid  element  [8]  the  ratio  of  inertia  forces  to  viscous  forces  is: 


and  the  ratio  of  inertia  forces  to  pressure  forces  is: 


V2  M2  YP/p  M2  yp 
Ap/p  ~  Ap  I ^  ~  Ap 


Therefore,  dividing  the  Reynolds  number  by  the  above  ratio  gives  the  ratio  of  pressure  forces  to  viscous  forces, 
R^Ap  .  It  should  be  noted  that  this  parameter  can  be  written  in  an  incompressible  form,  namely,  Re  t  where  the 
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pressure  coefficient  is  defined  in  the  usual  way  in  terms  of  the  reference  change  in  pressure  and  the  freestream 
dynamic  pressure. 

It  is  appropriate  to  investigate  the  data  used  to  develop  previous  correlations  and  understand  how  they  may  be 
predicted  from  correlations  based  on  the  new  parameter  for  laminar,  compressive  interactions. 

In  Figure  1  the  various  lengths  are  shown  for  a  two-dimensional  compression  ramp. 


Figure  1:  Nomenclature  for  Elliptic  Regions  Caused  by  2D  Compression  Ramp 

The  correlations  in  the  present  work  predict  the  ratio  of  the  upstream  influence  distance,  lu  ,  and  elliptic  distance 
(which  is  the  sum  of  the  upstream  and  downstream  distances),  l,  to  the  boundary-layer  thickness  SL.  The  distance  L 
is  defined  as  the  length  measured  from  the  leading  edge  to  the  source  of  the  disturbance.  In  the  case  of  an  impinging 
shock,  L  is  the  location  of  the  impingement  point.  Correlations  and  data  for  2D  laminar  compression  ramps  and 
shock  impingement  studies  were  published  previously  [9].  Additionally,  data  for  laminar  expansions  [10],  turbulent 
compressions  [11],  and  expansions  [12]  are  also  available  and  will  be  the  subject  of  future  work.  In  those  previous 
studies,  the  forms  of  the  correlations  were  developed  using  a  statistical  approach  combined  with  cross-plotting 
techniques  that  assumed  three  independent  parameters  for  adiabatic  flows:  Mach  number,  M,  Reynolds  number,  Re 
and  the  relative  change  in  pressure  in  the  inviscid  flow, (A p)/p_°°  .  The  present  study  treats  the  new  similarity 
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parameter  as  the  only  parameter  in  determining  streamwise  length  scales  for  elliptic  interactions  that  are  scaled  by 


the  undisturbed  boundary  layer  thickness. 


Figure  2  contains  the  data  from  previous  correlation  studies  for  the  upstream  influence  length. 


laminar  interactions  for  both  compression  ramps  and  shock  impingements[9]. 


Figure  2:  Correlation  Results  for  Upstream  Influence  due  to  2D  Compression  Ramps  and  Shock 

Impingement  Flowfields 


The  numerical  data  for  ramps  and  shock  impingements  alone  were  used  for  determining  the  relation  for  the  upstream 

1  f  Rg  Ap  \  2/5 

influence  in  terms  of  undisturbed  boundary  layer  thickness,  lu/SL  =  -I — ^ - )  .  The  experimental  data  for 

adiabatic  wall  conditions  was  taken  for  ramps  [13-16]  and  shock  impingements  [17-19]  from  previous  work.  The 

results  in  Figure  2  show  good  agreement  with  the  data.  In  this  figure,  the  correlation  value  for  the  upstream 

influence  is  defined  to  be  zero  when  the  change  in  pressure  is  zero.  Previous  correlations  were  developed  for  ramps 

and  shock  impingements  separately  for  Mach  number  and  Reynolds  number  ranges  of  2.0  to  6.06,  and  1  x  104to  8  x 

105,  respectively  [4].  Flowever  the  present  new  correlation  combines  both  ramps  and  shock  impingements  into  one 
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correlation.  Therefore,  the  previous  correlations  are  slightly  more  accurate  than  the  present  correlation  for  the 
upstream  influence  length.  Although  the  present  correlation  for  the  upstream  influence  length  is  not  as  accurate  as 
previously  developed,  the  results  for  the  total  length  or  elliptic  length  are  much  improved  as  shown  in  Figure  3. 


Figure  3:  Results  for  Elliptic  Region  Length 


Figure  3  has  the  data  used  to  develop  previous  correlations  [9]  for  the  length  of  the  elliptic 

(R.6  Ap  \  0.46 

— 2 - )  ,  plotted  as  a  function  of  the  elliptic  parameter.  As  was  done  for  the  results  in 

y  MqqPoq/ 

Figure  2,  the  correlation  developed  in  Figure  3  was  developed  from  only  numerical  data  for  both  compression  ramps 
and  shock  impingements.  There  were  no  experimental  data  from  previous  references  available  for  the  total  elliptic 
length  as  discussed  in  Ref.  [9].  The  present  correlation  for  total  elliptic  length  fits  the  data  very  well  and  enables 
one  to  determine  a  priori  the  total  length  of  an  elliptic  interaction  for  2D,  laminar  interactions  over  compression 
ramps  and  shock  impingement  flowfields. 
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III.  Conclusion 


New  correlations  have  been  developed  for  2D,  laminar  adiabatic  wall  flows  over  compression  ramps  and  flows 
with  shock  impingements.  These  correlations  are  derived  from  existing  numerical  data  and  incorporate  the  use  of  a 
new  similarity  parameter  for  the  Navier-Stokes  equations  that  represents  the  ratio  of  pressure  forces  to  viscous 
forces.  The  new  correlations  are  in  excellent  agreement  with  previously  published  data.  Development  of 
experimentally-derived  correlations  for  transitional  flows  over  compression  ramps  is  currently  underway. 
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